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Abstract 



Underwater communication is a challenging topic due to its singular 
channel characteristics. Most protocols used in terrestrial wireless com- 
munication can not be directly applied in the underwater world. In this 
paper, we focus on the issue of energy efficient transmission in under- 
water sensor networks (UWSNs) and analyze this problem in a rigorous 
and theoretical way. We formalize an optimization problem which aims 
i-^i to minimize energy consumption and simultaneously accounts for other 

performance metrics such as the data reliability and the communication 
r/3 delay. With the help of Karush-Kuhn- Tucker conditions (KKT condi- 

tions), we derive a simple and explicit, but nevertheless accurate, approx- 
imate solution under reasonable assumptions. This approximate solution 
provides theoretical guidelines for designing durable and reliable UWSNs. 
Our result also shows that reliability and communication delay are crucial 
factors to the energy consumption for transmission. 

Energy efficient transmission, underwater sensor networks, 
iS" KKT conditions. 

q 

O 1 Introduction 

< ^> Understanding the key mechanisms of the oceans is crucial to the knowledge 

of our Earth's climate and atmosphere. Over the past few years, there has 
T been a relentless effort to investigate this underwater world which remains little 

rS explored. This infatuation is highly motivated by various commercial, scien- 

tific, and military applications. The underwater sensor network (UWSN), a 
promising technique for the human being to explore the oceans, has a myriad 
of applications such as mineral exploitation, environmental monitoring, disaster 
prevention, military surveillance and coastline protection. 

Communication in underwater is a very challenging topic which is featured 
by large propagation delay, limited available bandwidth, high error probability 
and severe energy limitation. These features make the design of transmission 
mechanisms rather awkward. Among these features, energy efficiency is critical 
for a UWSN to maximize its utility. 

UWSNs are typically characterized by severely energy-limited nodes. For 
one thing, the nodes have small batteries and therefore cannot afford much 
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energy to complete tasks. For another, their power supplies cannot typically be 
replenished; that is, once exhausted, the node is discarded. Thus, efficient use 
of energy is of great importance to long-sustained and well-operated UWSNs. 
In this paper, we focus on the issue of energy efficient transmission. 

The issue of energy efficiency is widely considered in the literature. Various 
data transmission mechanisms [1-6] are proposed to achieve the objective of 
energy efficiency. These efforts are of great value to reduce energy consumption 
in UWSNs. However these mechanisms are designed more or less based on 
heuristics, and a theoretical analysis of energy consumption in UWSNs is little 
addressed. 

In this paper, we analyze energy efficient transmission in a rigorous and 
theoretical way. We formalize an optimization problem which aims to mini- 
mize energy consumption but also comprehensively considers other performance 
metrics such as data reliability and communication delay Using this problem, 
we study the effect of four key parameters (node distance, communication fre- 
quency, packet length and SNR) on the performance of UWSNs. Utilizing a 
two-step approach, we simplify this problem, and derive an explicit approximate 
solution to the simplified problem according to Karush-Kuhn- Tucker conditions 
(KKT conditions) [7]. We show that the approximate solution is accurate with 
regard to the simplified problem, and it also achieves near optimum of the 
original optimization problem. This approximate solution provides theoretical 
guidelines for designing durable and reliable UWSNs. Our result also shows 
that reliability and communication delay are crucial factors to the energy con- 
sumption for transmission. 

The rest of the paper is organized as follows. Section [2] introduces the system 
model of our paper. Section [3] formalizes the optimization problem and reduces 
it to a simplified version. Section [4] solves the simplified problem with the help of 
KKT conditions. Numerical result is presented in Section|5] Finally conclusions 
are drawn in Section [6] 

2 Preliminaries 

2.1 Scenario 

In this paper, we consider the scenario as is depicted in Fig. [l] There is a 
source node and a destination node. And the distance between the source and 
the destination is d. The source node keeps sending a package repeatedly until 
the destination node receives it correctly. This scenario, simple though it is, can 
be regarded as a basic component of multi-hop scenarios. Moreover, by studying 
it thoroughly, we can get valuable insights into more sophisticated schemes such 
as Forward-Error-Correcting based data transport protocols. 
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2.2 Energy Model 

The passive sonar equation [8] indicates that the signal to noise ratio (SNR) per 
bit 7b of an emitted underwater signal at the receiver 

7b = SL-TL- NL + DI (1) 

where SL is the source level, TL the transmission loss, NL the noise level and 
DI the directivity index. The unit of all quantities is dB. Since omnidirectional 
hydrophones are used, DI is 0. 

The transmission loss TL over a distance d (in m) for a signal of frequency 
/ (in kHz) is 

TL= 101gd+a(/) x dx 1CT 3 (2) 

where a(f) is the absorption coefficient in dB/km. It follows from the Thorp's 
formula [9] that 



«(/) = + jTT^n + 2J5 x 10 ~ 4 f 2 + °- 003 ( 3 ) 



0.11 J 2 44f 2 
1 + / 2 + 4100 + / 2 

The noise level NL is affected by four sources: turbulence, shipping, waves, 
and thermal noise. For the sake of convenience, we simplify NL with a practical 
approximation [10] 

7VL = 50-181g/ (4) 
Combining Q, (^J, ^ and Q, we obtain 

SL = 50 + 7b + 10 lgd +a(/) X dx 10~ 3 - 18 lg/ (5) 

On the other hand, / t (W/m 2 ), the transmitted signal intensity at lm from 
the source, is expressed as follows 

It = 10 SL/1 ° x 0.67 x 10~ 18 (6) 

Thus, the transmitter power P t , which achieves intensity I t at a distance of lm 
from the source in the direction of the receiver, is 

Pt = 2nxHxI t (7) 

where P t is in W, and H is the water depth in m. 
Utilizing 

we can express 7b with P t , d, and / 
7b =10 x (lg P t - lg(2irH x 0.67 x 10~ 18 ) - lg d) 
-a(f) dx 10- 3 + 181g/ 
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2.3 Channel Model 



In our model, we employ binary phase shift keying (BPSK), one of the most 
commonly used modulation schemes. In BPSK, the SNR per symbol, j s , equals 
the SNR per bit, 7b, that is, 7 S = Y&. 

The multipath effect is distinct in the underwater channel, so a Raylcigh 
fading channel is an appropriate model for the underwater channel (both in 
shallow and deep water) . The average BER of BPSK in an underwater Rayleigh 
fading channel is [11] 

P 6 (7>H(1-^) (9) 
where % = lCT^ 10 = 1(W 10 . 



3 Problem Formulation 

We concentrate on the energy consumption and its relations to the four key 
parameters (node distance, communication frequency, packet length, and trans- 
mitter power) in UWSNs. We also take into account other performance metrics 
such as reliability and communication delay. 

We optimize the energy consumption through the metric E , average energy 
consumption for successfully transmitting one bit payload. E can be written 
as follows 

E b = (io) 

lace ' l^p 

where E per is the energy consumption per transmission attempt, P acc is the 
packet acceptance ratio, and L p is the length of the payload for one packet. 
Meanwhile, some constraints should be satisfied. First, the packet acceptance 
ratio Pace should be no less than a threshold P acc o to maintain the reliability 
In addition, a packet must contain payload, and the communication frequency 
/ and the transmitter power P t must be nonnegative. 
Thus, we come to an optimization problem 



t-^ Ep er 
mm h = 

Pt,L,f Pace ' Lp 

Pt>0 



subject to < 



L> fl+T 

-Pace ^ PaccQ 

/>o 



(ii) 



where L is the length of a packet, /i the length of the header, and r the length 
of the trailer of a packet. Thus L p = L — (/i + r) . The optimization problem 



(10) aims to minimize the energy consumption and maintains the reliability at 
the same time. Moreover, the constraint of P aC c is also responsible for the com- 
munication delay. The link delay takes a dominating part in the total time cost 
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for one transmission (we will discuss it in Section [5| . For a given distance d, 
the link delay between the source node and the destination node of one trans- 
mission attempt is fixed. Thus the total time cost for successfully transmitting 
one packet, which is mainly composed of the link delay, largely depends on the 
packet acceptance ratio P acc . 

Now we develop an expression of E using P tl L and /. We have known that 

L p = L-fi-T (12) 

Pace can be expressed as 

/ / 1076/10 \ L 

Pace = (1 - PM) L = 0.5 + 0.54 / - ; — m (13) 



1 + 1076/10 

where 7b can be calculated using (Jsj in terms of P t , d and /. 

Furthermore, we also assume that if the bandwidth is / kHz, the available 
bit rate is / kb/s. In addition, for commercial hydrophones, the received energy 
of each packet is typically around one fifth of the transmitted energy P t [12]. 
Besides, a sensor node also consumes some extra power each time it transmits 
or receives a packet, such as energy used by MPU/DPS, A/D or D/A converter, 
filter and so on. This energy consumption should also be considered and we 
denote it by P c . So the total energy consumption for one transmission attempt 
is 

^ = ToW x( ^ t + Pc) (U) 

where L is the packet length. 



3.1 A Simplification of (11) 



From (n2b-(n4|, we find that in the optimization problem (111, the terms which 



involve / are rather complicated. This incurs great difficulties in solving (111 



To reduce complexity, we adopt a two-step approach to analyze (10 1. We first 



find optimal communication frequency f*(d) which minimizes the transmission 



power P t with respect to d, and reduce ( 10 1 to an optimization problem of two 
variables (P t and L). Second, we solve the reduced optimization problem with 
the help of KKT conditions. 

From ^ and Q, we know that SL directly determines the transmission 
power P t and that minimizing P t is essentially minimizing SL. So we differen- 
tiate ([5]) with respect to / and obtain 



dSL _/2.2 x 10~ 4 360.8 
~df-{ (l + /2)2 + (4100 + f 2 ) 

+ 5.5x10-)/^ — 



(15) 
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The node distance d in UWSNs is generally within the range of 0.1km and 
100km [13]. Moreover, it can be shown that for 0.1km < d < 100km 



,. dSL , ,. 6SL 

lim — - — = — oo, and lim — - — = +oo. 
f^0+ c)J f^+oo df 



(16) 



Thus it can be reasoned that Vd G [0.1,100], there exists an optimal fre- 
quency /*(d) which minimizes SL and hence Pt, and that ^§j ± \j—f (<j) = 0- So 
we can find /*(d) by letting = an d numerically solving the equation. The 
relation between /*(d) and d is illustrated in Fig. [2j 
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Figure 2: f*(d) vs. d. 



Substituting /* (d) for d, we can reduce (111 to a problem with two variables 
which has a relatively plain form 



min El 

Pt.L 



C 2 L(5P C + 6P t ) 



CxP t 



(17) 



' :/ - /' r) (0.5 + 0.5^ Cii+( . ,., 
P t > 

subject to ^ L>/i + r+ l 

accO 

where C* = 4.10 x lO" 13 x 10 10 ~ idxa (Md))dH, C x = /*(d) 18 and C 2 = 5000 ^ (d) . 
Note that Co, Ci and C2 vary with d, and they can be regarded as constants 
for a fixed d. 

For the convenience of further mathematical manipulation, we use the nat- 
ural logarithm of E' b as the objective function instead of E' b , which does not 
change the optimal solution. Thus, we obtain 

min In E' b = In C 2 + InL + ln(5P c + 6P t )- ln{L — fi — r) 



Lln(0.5 + 0.5 



CiP t 



C + C x p 



L> /i + t + 
) subject to { P acc > P acc0 
Pt>0 



(18) 



After exploring qualitative properties of (18 1, we will solve this problem using 
KKT conditions. In our model, we set H = 10m, P c = 1//W, and /i = t = 
16bits. 
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4 Finding Minimal Energy Consumption 



4.1 An Analysis of the Constraints 



The feasible region of (U8| is determined by the three inequalities in (18 1. How- 



ever, we will show that one of them is redundant, so they can be further sim- 
plified. 

Note that 



Pa. 



0.5 + 0.5 



CiP t 



(19) 



C + CiP t , 

We let 0.5 < P aC cO < 1- which is a reasonable assumption, so the inequality 
> P acc o can be rewritten as 



P 



P t > ^- 1 + 2 ^o) 2 >Q 



(20) 



^ClPaccoi 1 ~ Paced) 

Thus, we can remove Pt > from the set of constraints without changing the 
feasible region of the original optimization problem. 



According to (|20|) , we obtain an equivalent set of constraints 

C (-l + 2P a t c0 ) 2 



P 



> 



(21) 



4C 'l- P ac C o( 1 _ Pacco) 

L-(ji + t + 1)>0 
We denote the first inequality by hi(P t ,L) and the second inequality by 



li2(Pt,L). Then the optimization problem (18 1 is transformed to an equivalent 
one which has two inequality constraints 



f/n(P t ,L)>0 
min In Si, subject to < 
PfLeR \h 2 (P t ,L)>0 



(22) 



An illustration of the feasible region of ( 22 ) can be found in Fig. [3] (the hatched 
area). 



4.2 Qualitative Properties of (22) 



Before solving this problem, we first explore some of its qualitative properties. 



First of all, we explain some terminologies with respect to ( 22 1 . The feasible 
region is the set of points which satisfy all the constraints. We say that an 
inequality constraint hi(x) > is effective at a point x = (Pt, L) if the constraint 
holds with equality at x, that is, we have hi(x) = 0. E denotes the set of effective 
constraints at x. Obviously, E C {/ii,/i2}. If D(E) — {Vhi(x),hi e E} is 
linearly independent, we say that the constraint qualification is met at x. Note 
that V is the gradient operator. And if E = 0, we say that the constraint 
qualification also holds. 

The following two lemmas are of great value for finding the optimal solution 
to p2l. 
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Lemma 1. The optimization problem (22) satisfies the constraint qualification 
at any x — (Pt,L) within the feasible region. 



Proof. Note that 



V/ii(x) = (1, 



4L 2 Ci (- 



1 i pL \ 2 

1 ' r accO 1 



[Vh 2 (x) = (0,1) 5 



and that D(E) C {V/iiO), V/i 2 (x)}. If P>(P) ^ 0, it is obvious that D(£) 
is linearly independent, so the constraint qualification holds. If D(E) = 0, 
since the constraint qualification pertains only to effective constraints, it holds 
vacuously in this case also. □ 



Lemma 2. A global minimum exists to the optimization problem (22) 



Proof. In the expression of In E' b , we have In L > ln(i — fi — r) , VL > fi + r + 1, 
and ln(0.5 + 0-5\/ c C +cp ) < 0, VP* > 0, so we have 



In > ln(5P c + 6P t ) + In C* 2 



(23) 



For a given (Pto, Lq) in the feasible region, there exists an Mp t such that Mp t > 
P w and that VP t > M Pt , ln(5P c + 6P t ) + In C* 2 > In E' b (P m , L ). Using we 
know that In ^(P t) L) > ln(5P c + 6P t ) + In C 2 > In ^(P t0 , L ), VP t > M Pt . 



Thus, the feasible region of (22 1 can be divided into two areas by Mp t , as 
is shown in Fig. [3] Note that (P t Q,L ) lies within Area I. Moreover, accord- 
ing to the previous analysis, we know that V(P t ,i) which lies within Area II, 
\aE' b {P u L) >\nE' h {P tQ ,L ). 

On the other hand, Area I with its border is bounded and closed, so it 
is compact and the Weierstrass Theorem guarantees that ln_E£ has a mini- 
mum (P*,L*) in Area I with its border. Given the fact (P t o,Lo) lies within 
Area I, we can see that In P£(P t0 , L ) > \nE' b {P* ,L*). Since \nE' b {P t ,L) > 
lnE' b (Pto, Lq), V(Pt, L) which lies within Area II, we conclude convincingly that 
(P t * , L* ) is a global minimum to the optimization problem (22 1. □ 



After proving the above two lemmas, we introduce the following theorem, 
which claims the validity of our method. 

Theorem 1. Suppose the following conditions hold: 



1. A global minimum x* — (P t *,£*) exists to the optimization problem (22) 

2. The constraint qualification is met at x* . 



Then, the global minimum x* is a solution to KKT conditions of (22). 

Proof. Due to a lack of space, please refer to [7]. □ 
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I 



Pt = 



C (-l+2P^ c0 ) 2 
-lC 1 F a ^ c0 (l-P o ^ cQ ) 







'P„,.L„) 




/ I 





Figure 3: The feasible region (the hatched area) of (22 1 



4.3 Finding the optimal solution using KKT conditions 



We have shown that the optimization problem ( 22 1 has a global minimum and 
that the constraint qualification is met at the minimum. According to Theorem 
[l] we can find the optimal solution by solving KKT conditions. The KKT 
conditions [7] of (22 1 are listed as follows. 



d\nE' b 

dP t 
dlnE' b 
dL 



Ai 



dhi 

m 



x 



dh 2 
'dP 







dhx dh 2 _ 
Xl W X2 ^L-° 
hi(P t ,L)>0 

h 2 (Pt,L)>0 
\ 1 h 1 (P t ,L) = 
X 2 h 2 {P t> L) = 
Ai,A 2 > 



Equations (24al-(24gl are solved through a case by case analysis. 



4.3.1 Case 1, Ai ^ 0,A 2 ^ 



Pt = 



L = fi 



1 Op 1 + M + t 

1 AIr atxQ 



4c 1 (i-p a i +r 



(24a) 

(24b) 

(24c) 
(24d) 
(24b) 
(24f) 
(24g) 



In this case, ( |24e) and ( |24f) indicate that hi(P t ,L) = 0,h 2 (P t ,L) = 0. Thus, 



(25) 
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The objective function 



In El = lnC 2 + ln(5P c - 



3C (1 - 2Yf 
2C 1 {-1 + Y)Y 



(26) 



ln(Y 



— 1 — fl — T 



) + ln(l + /i + r) 



where Y = P££* 



In this case, (P t ,L) obtained from (251 is already a feasible solution to (22l. 



4.3.2 Case 2, Ai ^ 0, A 2 







In this case, (24el requires that 

P t = 

that is 



C (-l + 2P a ^) 2 

4ClP to(l --Palo)' 



(27) 



PaccO — -Pace — 0.5 + 0.5 



Co + C Y P t 



(28) 



Combining (27 1, (24a I and (24b), we can obtain an equation of X (X 



CiPt+c ) whose intricate appearance makes its solution a formidable task. 
So we follow an intuitive and pellucid way to tackle this problem. Observing 
that PaccO « 1 and 4- « 0, we approximate PaccO l to 



PaccO ^ 



1 



ln(PaccO) 



Substitute (|29| into (|27j, we get 

Co /ln(Pacco) 



p 



4Ci 



In ( PaccO ) 



(29) 



(30) 



Here it is more convenient to take E' b as the objective function. So we apply 
(29) and (28 1 to the expression of E' b in (17), and get 



C 2 P(5P c + 6P t ) 
L — fl — T 
1 



PaccO 

c 2 



A 



accO hp 



(— + 3C L P ) + P 



(31) 



where L r , 



PaccO \ 2Ci In P a , 

L-fi-r. In Q, 

A =lnP QCc0 ( M + r)(9C - 10dP c ) + 3C (ln P acc0 ) 2 
+ 3C (Ai + T) 2 

lnP acc0 (10CiP c - 9C ) - 6C (m + r) 



2C 1 lnP a 



ccO 
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It can be showed that lnP acc o < and A > 0, so when 



L = /! 



A 



(32) 



E' b , or ]nE' b , achieves its minimum. The corresponding optimal Pt can be com- 
puted using ( [30| . 

In this case, we attain an explicit approximation instead of an exact answer. 
In Section [5] we will show that this approximation receives satisfactory results. 
More specifically, given a wide range of d and a practical value of P aC cO, the 



(Pt,L) obtained from (30) and (32) can meet the threshold packet acceptance 
ratio PaccO and save considerable energy. 



4.3.3 Case 3, A : = 0, A 2 ^ 

In this case, L = /i + r + 1 and 
and simplifying it, we obtain 



9 In E' b 
dPt 



0. Applying L = /i + r + 1 to (24a I 



a + a x X + a 2 X 2 + a 3 X A = 



(33) 



where X = J „ „ and X £ (0, 1). The coefficients of (33) are 



Ca+CtPt 

ao = hC\P c L 

a 2 = (6C - bdP c )L 



a i 
a?, 



-12C - BdPcL 
-(6C -5CiP c )i 



We denote the left part of (|33| by T 3 (X). Note that T 3 (0+) = 5CiP c (l+/i+r) > 
and that T 3 (l") = -12C < 0. Thus ((33} has a root within (0, 1) which can 



be found by numerical solvers such as Bisection method. 
4.3.4 Case 4, Ai = 0, A 2 = 



Combining ( 24a I and ( 24b I , we obtain 



where 



X = 



fi(X) 



CjPt 
Co+dPt 



f 1 (X) + f 2 (X) 



1 24C X 

2 (-1 + X) (-6C X 2 + 5dP c (-1 + X 2 )) 

+ r)(ln 16 + {n + r)(ln(i±^)) 2 - 41n(l + X)) 
21n(¥) 



(34) 



and X e (0, 1). We denote the left part of (34 1 by T 4 (X). Note 



that T 4 (0+) < and that T 4 (l~) = +oo. Thus, (34 1 has a root within (0,1) 
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5 Numerical Result 



In this section, we assess the solutions to the four cases in Section [4j and obtain 
the conclusion that ( 30 ) and ( 32 ) constitute an accurate approximate solution 
to (22 1 given a considerably wide range of d and a practical value of P aC cO- We 



also show that our two-step approach to analyze (111 can effectively find near 



optimum of (111. 



5.1 Feasibility of Solutions to (22) 



In Section |4j we obtain four solutions (four cases) to (22 1. Here, we check 
whether these solutions satisfy two constraints (hi and h^) of (22 1. It can be 



shown that for each case, the constraint of L (/12) is met. Thus we investigate 
the packet acceptance ratio P acc of each solution. 

For Case 3 and Case 4, Ai = 0, so the solutions are independent of P aC cO 
and hi might not be effective in Case 3 and Case 4. Fig. [4] depicts the packet 
acceptance ratio of solutions to Case 3 and Case 4. We can see that in both 
cases, P acc decays fast as the node distance d grows. On the other hand, a 
typical threshold P aC cO is usually no less than 95%. So the solutions to these two 
cases are not feasible except for very small node distance. This indicates that 
reliability and communication delay are crucial factors to the energy efficient 
transmission in UWSNs. In other words, considering energy consumption alone 
incurs high packet loss and large communication delay. 




node distance d (in in) 



Figure 4: Packet acceptance ratio of solutions to Case 3 and Case 4. 



According to ( 25 1 , the solution to Case 1 is a feasible solution to ( 22 ) (Fig 



|5j). Although in Case 2, we derive an approximate solution instead of an exact 
solution , we find in Fig. [5] that P acc of the approximate solution converges to 
PaccO fast as d grows. In addition, for a higher P aC cO, the converge becomes 
faster. In other words, the approximate solution to Case 2 satisfies hi for a 
sufficiently large d and a sufficiently high P acc o- 



5.2 Optimality of the Approximate Solution 

According to the above analysis, we find that the solution to Case 1, and the 
approximate solution to Case 2 are feasible with respect to hi and h 2 - In Fig. 
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0.5 1 1.5 2 2.5 3 3.5 4 

node distance d (in m) x 10 4 



Figure 5: Packet acceptance ratio of solutions to Case 1 and Case 2 with respect 
to different P aC cO- 



[6j we see that the value of objective function corresponding to the solution to 
Case 1 lies in the upper part of the graph. Thus the approximate solution 
attains better performance than the solution to Case 1. We also numerically 
solve (p2|, and compare the numerical solution with the approximate solution 
(TableTiJ. We convincingly conclude that the approximate solution to Case 2 



is an accurate approximation of the solution to ( 22 ) 



\ d 


10000 


20000 


30000 


40000 


50000 


0.980 


0.012 


0.015 


0.015 


0.015 


0.015 


0.985 


0.009 


0.011 


0.011 


0.012 


0.012 


0.990 


0.007 


0.000 


0.000 


0.000 


0.000 


\ d 


60000 


70000 


80000 


90000 


100000 


0.980 


0.015 


0.015 


0.015 


0.015 


0.015 


0.985 


0.012 


0.012 


0.012 


0.012 


0.012 


0.990 


0.000 


0.000 


0.000 


0.000 


0.000 



Table 1: Relative error (%) of the approximate solution compared to the nu- 
merical solution to ( 22 ) 



In Section |4| w e adopt a two-step approach to analyze the original optimiza- 
tion problem (111. Now we numerically solve (1 1 1 h , and compare the numerical 



solution with the approximate solution in Fig. km We see that the approximate 



solution we derive can find near optimum of the original problem (111. 

5.3 Verification of the Assumption in Section [3] 

We verify the assumption made in Section [3] which states that the link delay 
takes a dominating part in the total time cost for one transmission. We denote 



13 



Figure 6: Comparison of the objective functions of different solutions. The 
circle marked data represents the solution to Case 1, the square marked data 
represents the approximate solution to Case 2, and the triangle marked data 
represents the numerical optimal solution to (111. Different colors correspond 
to different P acc o- 



the link delay by T\ and 



where d is the node distance and v is the propagation speed of acoustic signals 
in water (here we let v = 1500m/s). On the other hand, the sending or receiving 
time for a packet, T2, can be expressed as 



To = 



L 



1000/ 

where L is the optimal packet length in bit and / is the optimal frequency in 
kHz. From Fig. [7J we find that T\ ^> T%. Thus, the assumption is validated. 




Figure 7: |f vs. 
same. 



d. For different P aC cO, the shape of the curve is almost the 



5.4 Summary 



In this section, we demonstrate that (30 1 and (32 1 constitute an accurate ap- 
proximate solution to (22) given a wide range of parameters. We also show that 
( 30 ) and ( 32 1 can find near optimum of (111. 
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We vary d and P aC cO, calculate the corresponding optimum value of hiE' b 
using the explicit approximation (30 1 and (32 1, and obtain Fig. [8] 




Figure 8: The optimum of (22 1. 



From Fig. [6] and [8] we can see that the energy consumption grows as P aC cO 
and d become larger. This is consistent with intuition. When the node distance 
becomes larger, more energy is required keep a reliable transmission. On the 
other hand, more energy is necessary to meet higher P aC cO- 

From Subsection 



5.2[ we know that hi is an effective constraint for the 
optimal solution to ( 22 1 . This indicates that reliability and communication 



delay are crucial factors to the energy consumption for transmission. That is, 
concerning reliability and delay is a necessity when we intend to minimize energy 
consumption for transmission. 



6 Conclusion 

Following a rigorous and theoretical way, we address the issue of energy efficient 
transmission in UWSNs in this paper. Firstly we formalize an optimization 
problem which aims to minimize energy consumption but also comprehensively 
considers other performance metrics such as data reliability and communication 
delay. Using this problem, we study the effect of four key parameters (node 
distance, communication frequency, packet length and SNR) on the performance 
of UWSNs. Then we derive an explicit approximate solution to the optimization 
problem according to KKT conditions. This solution can provide theoretical 
guidelines for designing durable and reliable UWSNs and we show that the 
approximate solution is accurate with regard to a wide range of parameters. Our 
result also shows that reliability and communication delay are crucial factors to 
the energy consumption for transmission. 
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